A=[ 1 3 3 4 5];
B=[6 7 8 9 10];
C = num2cell([A;B],1)

mat1=[
    c2,-c1;
    -b2,b1];
vec1=[x3;y3];

cellfun(@(para)(mat1*para+vec1)',C)
